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Abstract 

wo: 

Starting from the study of the low-energy and high-energy behaviours of the QCD 
\ three-point functions (VAP), (VVP) and (AAP), several 0(j? 6 ) low-energy constants of 

^r) ' the chiral Lagrangian are evaluated within the framework of the lowest meson dominance 

(LMD) approximation to the large-Ac limit of QCD. In certain cases, values that differ 
substantially from estimates based on a resonance Lagrangian are obtained. It is pointed 
out that the differences arise through the fact that QCD short-distance constraints are in 
general not correctly taken into account in the approaches using resonance Lagrangians. 
We discuss the implications of our results for the C(p 6 ) counterterm contributions to the 
\& | vector form factor of the pion and to the decay n — > ef e 7, and for the pion-photon-photon 

transition form factor. 



"p^' 1 Introduction 

^| In the chiral limit, the lightest pseudoscalar states of the hadronic spectrum become the octet 
of massless Goldstone bosons resulting from the spontaneous breaking of the chiral SU(3)l x 
SU(3)r global symmetry of the QCD Lagrangian towards its diagonal subgroup SU(3)v of 
vector symmetries. This well-known fact [I], allows to describe the interactions of the light 
pseudoscalar mesons at low energies in terms of an effective Lagrangian 1, 1. The latter 
involves the pseudoscalar fields, described by a unitary matrix U(x), and transforming under a 
nonlinear representation of the chiral symmetry group, as well as the sources, v^x), a^x), of the 
light-quark vector and axial currents, and s(x), p(x), of the scalar and pseudoscalar densities 
of QCD ||]. Matrix elements of these currents between pseudoscalar states, or scattering 
amplitudes involving these light states only, can be computed in a systematic way in the low- 
energy theory as long as all momentum transfers p 2 are sufficiently small, p 2 <C h 2 H , where Ah ~ 
1 GeV is the typical scale at which the non-Goldstone bound states of QCD are formed. Since 
the (running) light quark masses m q (u) are also small as compared to this scale, m Ui d jS (A#) <C 
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Ah, the effective Lagrangian C eS can be organized as an expansion in powers of derivatives of 
the field U(x) and powers of the light quark masses, 

k,l 

with 

p \ k ( m q 



£ <-> ~ it) ■ (L1) 

The presently available studies on the structure of the low-energy effective Lagrangian involve, 
beyond the lowest order terms, the pieces £(4,0), £(2,1) an< ^ A°> 2 ) @' A 6 >°) anc ^ A 4 . 1 ) & 0' A 2 , 2 ) 
and £(0,3) [0, ||, 0, as well as £(0,4) 0- The structure of each £(k,i) is entirely fixed by the chiral 
symmetry properties of QCD ||, but involves coefficients, the so called low-energy constants, Q 
which are not determined from symmetry requirements alone. The predictive power of the 
effective theory therefore hinges to quite some extent on the knowledge of these low-energy 
constants. At 0(p 4 ), the values of most of the low-energy constants were extracted from data. 
The proliferation of low-energy constants at 0(p 6 ) makes such an approach unrealistic. 

On the other hand, it is a general property that these low-energy constants correspond to 
the coefficients of the Taylor expansion, with respect to the momenta, of some QCD correlation 
functions, once the singularities (poles and discontinuities) associated with the contributions 
of low-momentum pseudoscalar intermediate states have been subtracted. The characteristic 
feature of the Green's functions that are actually involved is that they are order parameters 
of the spontaneous breaking of chiral symmetry. Thus, they do not receive contributions from 
perturbative QCD at large momentum transfers, but rather exhibit a smooth behaviour at 
short distances. The low-energy constants are thus expected to be sensitive to the physics 
in the intermediate energy region, that is to the spectrum of mesonic resonances in the mass 
region around the hadronic scale Ah- This basic observation underlies, in some way or another, 
most attempts to estimate the values of the low-energy constants from resonance data (for an 
introduction to the vast bibliography on this subject, we refer the reader to the review articles 



m mm)- 

It has become customary to describe the effects of resonance states within a Lagrangian 
framework, by introducing, besides the Goldstone boson field U(x), additional local fields as- 



sociated with the meson states. While there exists a systematic way |13] to construct fields 
which have the appropriate transformation properties under chiral transformations and invari- 
ant Lagrangians, there is however no restriction from chiral symmetry as to the number of fields 
and the order of derivatives thereof involved in the terms which describe the interaction among 
resonances or between the resonances and the Goldstone bosons. In addition, the construction 



of Ref. [13] leaves open the question of the choice of the Lorentz group representation for the 



field describing a meson state of a given spin. This lack of restrictions has led to many proposals 

1 Some of the counterterms involving the external sources only and no pseudoscalar fields actually rather 
correspond to "high-energy constants", since they describe the (perturbative) short-distance ambiguities of 
some QCD correlators. At 0(j? 4 ), this concerns the constants Hi and H2 of Ref. Q. We exclude this type of 
counterterms from the present discussion. 
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concerning a Lagrangian description of interacting Goldstone bosons and mesons (the litera- 
ture on this subject can be traced back from several reviews [TF|, [TR ITB[ and from the articles 



17| |18[ |19j). It has been pointed out in Ref. |2(J that restrictions actually can be introduced, 



via the QCD short- distance properties of the relevant Green's functions or form factors, and 
by requiring that these properties are satisfied by the same objects constructed with the help 
of the Lagrangian involving resonances. This aspect has been taken into account in a rather 
systematic manner at the level of the 0(p 4 ) counterterms |[20|| . Although resonance estimates 
have been given for several 0(p 6 ) low-energy constants, the importance of implementing the 
appropriate QCD short-distance constraints has not always been stressed. The purpose of the 
present study is precisely to address this aspect. 

Although quite practical and useful, since it guarantees that some general properties (lo- 
cality, analyticity) of quantum field theories are correctly taken into account, a Lagrangian 
formalism is not absolutely necessary in order to estimate the contributions of the QCD reso- 
nance states to the low-energy constants. In this article, we shall in fact consider a different 
approach, working directly with the appropriate Green's functions. We shall however retain 
the general features and assumptions that underlie the Lagrangian approach, and that we 
briefly recall. First, one usually considers zero- width resonances, and takes into account the 
contributions from one-resonance states, which produce only poles in the corresponding corre- 
lation functions. There exists a well-defined framework where this kind of restrictions arises 
in a natural way [^T], namely the large- Nc limit of QCD E3|. On the other hand, working 
in this limit still requires to consider, in each channel, an infinite number of resonances, with 
masses and coupling constants adjusted such as to reproduce the QCD perturbative contin- 
uum at high momentum transfers. Of course, such a point of view is rather ambitious, since 
it amounts to solving QCD in the large- Nc limit, a reputedly difficult task. We shall adopt 
a more modest attitude, assuming that in each channel a few lowest-lying resonances give al- 
ready the main contribution (approaches involving an infinite number of zero-width resonances, 
with various additional assumptions about their spacing, decay constants, etc., can be found 



in Refs. p3| , pl| f25fl ). The number of resonances to be considered in each channel will be taken 
as the minimal (finite) number necessary to satisfy the requirements set by, say, the leading 
QCD short- distance constraints for the Green's functions under consideration, and possibly 
other constraints that one may wish to impose. This minimal hadronic ansatz (MHA) ap- 
proximation may be well justified in the case of Green's functions which are order parameters, 
free of perturbative contributions (for Green's functions which are not order parameters, the 
QCD continuum contribution has to be included as well, see Refs. |26|, P^j). In fact, in many 
cases, this minimal hadronic ansatz can be reduced to retaining, in each channel, a single res- 
onance. At the 0(p 4 ) level, this lowest meson dominance (LMD) approximation to large-A^ 
QCD has been tested successfully in several instances |L7|, [L8|, |20], |26fl . As we shall see in the 
examples treated in the present work, depending on the constraints one wishes to implement 
on the Green's functions under consideration, this simplest LMD approximation may however 
not always be sufficient. A second feature common to the Lagrangian and to the LMD or MHA 
approximations is the fact that the estimates of the low-energy constants do not reproduce 
their scale dependence. The latter, which comes from Goldstone boson loop contributions to 
the relevant Green's functions, is a next-to-leading order effect in the 1/Nc expansion, and 
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lies thus beyond the approximation considered. We shall adopt the usual point of view that 
the estimates furnished by this type of approach corresponds to the values of the low-energy 



constants at the typical scale, say A#, set by these resonance states [IS 



Approaches which do not rely on a resonance Lagrangian have been used before at the 0(p 6 ) 
level for two-point functions [^, |29|, |30|. In these studies, the relevant low-energy constants 
have often been expressed through (super convergent) dispersion integrals of the corresponding 
spectral densities, which were then evaluated using available data. It seems difficult to follow 
similar lines in the case of three-point or higher Green's functions. Not only are their analyticity 
properties far more complicated, but the corresponding spectral densities are in general not 
experimental observables. Studies of three-point functions similar to the lines we follow here 
can be found in Refs. [[H], [32|, [33|], although the discussion of their short-distance properties is 
less complete than the one presented below. 

In this article, we shall concentrate on a certain subset of 0(p 6 ) low-energy constants con- 
tributing to £(4,1) and £( 6i0 ), an d corresponding to the three-point functions (VAP), (VVP) 
and (AAP) (see the beginning of section 2 for the precise definitions). There are several reasons 
for this specific choice. First, these correlators have been considered before in the literature, 
so that quite some information concerning them is already available. Second, although the 
approach that we shall follow here is, in principle, applicable to other correlation functions as 
well, these Green's functions have a certain degree of simplicity, which makes them particularly 
valuable for illustrating our point of view. Finally, these Greens functions also play a role 
in the evaluation of some of the counterterms that arise in the calculation of electromagnetic 
contributions to the pseudoscalar masses [34] or of radiative corrections to semileptonic decays 
of the pseudoscalar mesons |35|, |36|] within an effective Lagrangian framework []37| . 



The remaining material of the present article is organized as follows. In section 2, we 
define the relevant QCD Green's function and study their long-distance properties in the chiral 
limit and at leading order in the 1/Nc expansion. In particular, we identify the low-energy 
constants related to these correlators. Section 3 is devoted to an extensive discussion of the 
leading short- distance properties of these Green's functions within the same framework. In 
section 4, we construct some simple ansatze, in terms of a finite number of narrow resonances, 
which correctly reproduce the short-distance constraints. These are used to determine the 
corresponding low-energy constants in section 5. We then compare our results with those 



obtained from a Lagrangian involving resonances |38| which has often been employed in the 
literature to estimate the low-energy constants at order p 6 (section 6), and point out that this 
resonance Lagrangian does not correctly incorporate the necessary short- distance properties 
(section 7). In section 8, we present several applications. Conclusions and additional discussions 
can be found in section 9. Appendix A contains some technical details relevant for the discussion 



in section 4 and Appendix B gives the expression of the resonance Lagrangian of Ref. [38 
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2 Long-distance properties from chiral symmetry 

We consider, in the three flavour chiral limit, the momentum space QCD three-point functions 
(IW) £(p, q) = [d 4 x [ d 4 ye^ + ^\0\T{V^x)Ai(y)P c (0)}\(}} , 



(U VVP )f u c ip,q) = Jd'xJ d 4 ye^ x ^(0\T{V;(x)V u b (y)P%0)}\0) , 

(U AAP )^(p,q) = Jd'xJ rfV (p " + ^ ) (0|T{A«( a :)4(y)P c (0)}|0), 
involving the octet vector and axial currents, 

Afa) = (^75yV0(z), 
as well as the octet pseudoscalar density 

P a (x) = (^ 7 5y^)(x). 
These three-point functions satisfy the following chiral Ward identities, 



(2.1) 



p^U VAP )^ c (p,q) 



W)o/ 



abc 



<h _ (p + q)u 

q 2 [p + q) 2 \ ) 



?"(n 



VAP, 



\abc 



(M)nf abc ( p + q ^ 



^(rw)"t c (p, g) = o , g"(iwC(p, g) = o , 



abc I 



p^U AAP )f v c (p, q) = 0, q^U AAP ) abc (p, q) = , 



/IV 



(2.2) 



where (ipip)o denotes the single flavour bilinear quark condensate in the chiral limit. The 
general solution of these Ward identities, taking into account the invariances of QCD under 
SU(3)v, parity and time reversal transformations (the latter being responsible for the absence 
of structures of the type d abc in the case of Hvap, or of structures of the type f abc in the cases 
of Hvvp an d Tl AAP ), read 

(p + 2q) ll q l/ r]^ 



(IWff(p,g) 



(U VVP )p{p,q) 
(iWfffeg) 



f \ m \q^ + q? ( P + qr 
+P^(p, q)F(p\ q\ (p + qf) + Q^(p, q)G(p 2 , q\ (p + g) 2 )| , 
e^p a q^d abc H v (p 2 ,q 2 ,(p + q) 2 ), 

e,ua0 a q^d abc n A (p 2 ,q 2 ,(p + q) 2 ). 
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(2.3) 



5 5 5 



and we 

0^,1^,2^,3 



Here, rj^ denotes the metric tensor in flat Minkowski space with signature (+ 
use the conventions £0123 = 1 for the totally antisymmetric tensor e^ upa and 75 = i7°7 1 7 2 7 3 . 
The transverse tensors P^ and Q^ u are defined by 

Ppvfa l) = %Vv ~ (p • q)v»v > Q^ip, q) = V 2 %1v + q 2 p^Pv - (p • q)p^Qu - p 2 q 2 Vtiu ■ 
Due to Bose-Einstein symmetry, the invariant functions Tiv,A have the property 

HvXp 2 , q 2 , (p + q) 2 ) = HvA?y> (p + <?) 2 ) • ( 2 - 4 ) 

The behaviour of these invariant functions at small momentum transfers is constrained by 
the presence of singularities arising from Goldstone boson intermediate states. Here we are 
interested in the limit where the number of colours Nc becomes infinite. In this limit, the 
contributions from one-particle intermediate states dominate, so that at low energies we only 
need to keep the corresponding Goldstone boson poles and the polynomial terms involving the 
counterterms. In the even intrinsic parity case, we use the basis of Ref. for the 0(p 6 ) 
counterterms, and we obtain 



r 



ChPT/2 2 



(p\q\(p + qY) 



F${p + q) 2 



L 9 + L 
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5 2 1 2 

+ (CVs — 2^ 88 — ^ 89 + 3Cgo)p + (C 78 — 2C 87 + -C 88 )q 



+ (C 78 + 4C 82 --C 88 )(p + g) 2 



+ 



iChPTlJl 2 



(p ,q\(p + q)) 



F$q 2 (p + qy 

+2(-C 88 + C 9Q )p 2 + (2C 78 - C 89 + C 90 )q 2 - 2C 90 (p + qf 



+ 



(2.5) 



where the ellipses stand for higher order contributions. For the two correlators of odd intrinsic 
parity, we use the counterterm Lagrangian of Ref. || , in terms of which we find 



H C v hPT {p 2 ,q 2 ,{p + qY) 



W)o 
F$(p + q) 2 



N, 



a 



in' 



+ (AA 2 -16A 3 )(p 2 + q 2 



+{-4A 2 + 8A 3 + 16A 4 )ip + q y 



+ 



n c A hPT ( P 2 ,q 2 ,(p + q) 2 ) 



F^P + q) 2 



N, 



c 



247T 2 



+ (4A n + 4A 23 -16A 24 )(p 2 + q 2 



-UAu - I6A17 - 4A 23 + %A 2A + 16A 25 )(p + q)' 



+ 



(2.6) 



We have thus identified the set of low-energy constants that describe the long-distance behaviour 
of the (VAP), (VVP) and (AAP) three-point functions. 
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3 Short-distance analysis 



We next study the properties of the (VAP), (VVP) and (AAP) correlators at short distances. 
These will be conditioned by the fact that the three Green's functions under consideration 
are order parameters of chiral symmetry. Therefore, they vanish to all orders in perturbative 
QCD in the chiral limit, so that their behaviour at short distances is smoother than expected 
from naive power counting arguments. Two limits are of interest. In the first case, the two 
momenta become simultaneously large, which in position space describes the situation where 
the space-time arguments of the three operators tend towards the same point at the same rate. 
Our analysis will be restricted to the leading terms R, and the expressions below hold up to 
corrections of order 0(a s ). We obtain 



lim _(IW) f u c (Xp, Xq) 



A— >oo 



/ 



abc 



2„2 



P (P + 2 <?W - VfwP Q 



lim (n VVP ) a ;c(\p,\q) 

A— >oo ^ 



p 2 q 2 (p + q) 2 
\[p 2 -q 2 -{p + qf}P„-Q„} - 

q 2 +p 2 + (p + q)' 
p 2 q 2 (p + q) 2 



mbc 



d^e^^q 



hm(IW)J c (Ap,Ag) 

A— >oo ^ 



2A 2 

(Wo ,abc c « -9 q 2 +p 2 -{p + q) 2 
a ^iiuafiP q 



2X 2 



p 2 q 2 (p + q) 2 



+ 01* 



+ o - 



One concludes that 



lim TUXpf^Xqf^Xp + Xqf 

A— >oo 



lim G((Xp) 2 ,(Xq) 2 ,(Xp + Xqf 

A— >oo 



2 ,,2 (j) + q y 



r 



1 p 

2A?W>0- p 2 q 2 {p + q) 2 



A 6 p 2 q 2 (p + q) 2 



O 



A 8 



and 



hm H^((Ap) 2 ,(Ag) 2 ,(Ap + Ag) i 

A— >oo 



2A< 



T p 2 + q 2 ±(p + q) 2 



(3.1) 

(3.2) 
(3.3) 

(3.4) 



Notice that since the (V^o condensate and the pseudoscalar density P a (x) have the same 
anomalous dimensions, the leading short-distance behaviour exhibited in these expressions is 
canonical, the corresponding Wilson coefficients have no anomalous dimensions. 

The second situation of interest corresponds to the case where the relative distance between 
only two of the three operators involved becomes small. It so happens that the corresponding 
behaviours in momentum space involve, apart from the correlator (AP) which, in the chiral 
limit, is saturated by the single-pion intermediate state, 

d*xe ip - x (0\T{A a u (x)P b (0)}\0) = 5 ab (^) ^ , 

2 In the case of the (VAP) correlator, the subleading term in the short-distance expansion, involving the 
mixed condensate, can be found in Ref. |33[ . 
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the two-point function (VT) of the vector current and the antisymmetric tensor density, 

5 ab (U VT ), pa (p) = |d 4 xe^(0|T{\/;(x)(^a pCT ^)(0)}|0), 

with (7 pu = |[7p,7 CT ] (the similar correlator between the axial current and the tensor density 
vanishes as a consequence of invariance under charge conjugation). Conservation of the vector 
current and invariance under parity then give 

(TIvt)w*(p) = (PpVtia ~ PaVw) H-vt(p 2 ) ■ 
The leading short- distance behaviour of this two-point function reads (see also |39|] ) 

&nvT( (^) = 4^>£ + ^). (3,) 

For the (VAP) correlator, we then find 

Km (IW) f;(X P , q - Xp) = - I f abc (^) P^+P^~(P-g)^ + o (i) , (3.6) 

\abc(\^ „\ ^ eabc /J~,„l,\ PflQv , /n { 



and 



hm (n VAP ) J c (Ap, g) = -f^(W)o^ + ^ (jj) , (3.7) 



Km (IW) £ (p, Ag) = 1 r bc p » g " n yT (p 2 ) + O ) . (3.8) 



For the (VVP) correlator we obtain the results 

hm (U VVP )^(Xp, q-Xp) = - 1 rf afcc (^)o £fivpaP^Q cr + O (± ) , (3.9) 



A lim(n yy p)J c (Ap,g) = ^ abc e^ p(T ^n yT (g 2 ) + O(l) , (3.10) 
and for the (AAP) correlator, 



Urn (IW) f;(Xp, q-Xp) = - 1 <f bc (^>0 + O ) , (3.ir 



A lim(n AAP )J c (Ap,g) = Of^J . (3.12) 
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In terms of the invariant functions T and Q, the constraint 

(Wo " 



yields 



lim H{Xp)\{q-\p)\q^ 

A— >oo 



lim £((Ap) 2 ,(g-Ap) 2 ,g 2 ^ 

A— >oo 



X 2 p 2 

(Wo 
(AV) 2 



£ (( V) + 



A 



together with 



F {l \q 2 )-G (l \q 2 )+G m {q 2 ) 
Finally, the following properties 



lim T{{\p)\q\{q + \py 

A— >oo 



q- 
2 



lim £((Ap) 2 ,g 2 ,(g + A^ 

A— >oo 



o 
o 



1 

A3 
1 

A^ 1 



and 



lim -F(p 2 ,(Ag) 2 ,(p + Ag) 2 ) 



1 1 



A— >oo 



IlvAp 2 ) + o 



A 2 q 2 VA 3 



lim G(p 2 ,(\q) 2 ,(p + \q)' 

A— >oo 



(9 



A- 



1 

Y 2 



(3.13) 



(3.14) 



(3.15) 



(3.16) 



must also be satisfied. For the invariant functions Tiy and 'Ha we obtain the constraints 

1 ,- , , 1 



Urn H^((Ap) 2 ,(g-Ap) 2 ,g 2 ^ 

A^oo 



(3.17) 



lim n v ((XpY,q 2 ,(q + Xpy 

A— >oo 

lim n A ((\p) 2 ,q 2 ,(q + \pf 

A— >oo 



4 The intermediate energy region 



(3.18) 
(3.19) 



In this section, we shall construct representations of the invariant functions which describe the 
correlators (VAP), (VVP), (AAP) and (VT) in the intermediate energy region, dominated by 
the resonances, and which reproduce the short- distance constraints studied in the preceding 
section. Finding the general structure of the invariant functions J 7 , Q and 7~Cv,a is of course far 
beyond our present possibilities. As discussed in the introduction, we shall therefore work in 
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the framework of the large- Nq approximation to QCD and assume, in addition, that already a 
finite number of resonances will give a satisfactory description of these correlators. 

We first consider the case where only a single resonance is retained in each channel, assuming 
that in the pseudoscalar one only the massless Goldstone bosons need to be kept. The corre- 
sponding lowest meson dominance (LMD) ansatze for the invariant functions are constructed 
such that the resulting expressions agree with the short-distance constraints ( ft.2|) - Q3.4J) . This 
is rather straightforward, and in the case of the (VAP) correlator, the result reads 



F LMD (p 2 ,q 2 ,(p + q) 2 ) 



(ipip)o p 2 - q 2 - (p + q) 2 + 2a 

2 (p2_ M 2 )(g 2_ M 2 )(p + g) 



n 2 ^P + q) 2 ) = (#)o ^ M , )( ;; + M ^ +g)2 . (4.1] 



The constants a and b in Eq. ( |4.1| ) can be determined as follows. As shown in Ref. ||33|| , one can 
relate the (VA\tt) vertex function T V a (for the definitions and a discussion of some properties 
of the vertex functions associated to the three point functions under study, we refer the reader 
to Appendix A) to the two-point correlator (VV — AA) via a low-energy theorem. From the 



two Weinberg sum rules [HDJ one obtains in this way the relation 

a - b = -[Ml + M\) . (4.2) 

The constant b can be fixed by requiring that the vector form factor of the pion Fy(q 2 ), defined 
by 

W)IK b (o)K(p)> = *f abc (p' + p)»Fv(q 2 ), q = p-p', (4.3) 



satisfies an unsubtracted dispersion relation |2(| (there are theoretical arguments ^T, f|2| in 
favour of a 1/q 2 fall-off of this form factor for large momentum transfer) . Since 



Fy(q 2 ) = 1 + -i— lim lim (q -p) 2 p 2 g(q\ (q-p) 2 ,p 



we obtain with the function Q LMD from Eq. ( |4.1|) the result 

b q 2 
2M\ q 2 - Ml 

b = 2M\ . (4.5) 



F v (q ) = 1 " 77772 77 7F7 » ( 4 -4) 



and thus 

Combining Eqs. (^2) and (|4.5| ) then gives 



M 2 A -M^. (4.6) 



We note that the argument given in Ref. [33] that the same result for the constant b can also 



be obtained by enforcing the correct short- distance behaviour of the (V^PI^) vertex function 
Typ is not correct. In Appendix [A] we sketch the derivation of the operator product expansion 
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of this vertex function, with the result given in Eq. ( |A.6| ). The LMD ansatz will reproduce 
the subleading term in the OPE provided b = 4M|/3, which implies with Eq. ( |4.2[ ) a = 
—My + M^/3. Thus, with the LMD ansatz for the invariant function Q from Eq. ( |4.1|) it is 
not possible to reproduce at the same time the requirements from the asymptotic behaviour 
of the form factor Fy(q 2 ) and from the subleading terms in the OPE for the (VP\n) vertex 
function. The subleading term of the OPE of (VA\n) is correctly reproduced if Eq. ( |4.2|) holds. 
In what follows, we shall understand that the LMD approximation for the (VAP) correlator 
corresponds to the ansatz (|4.1| ) together with the choice b = 2M\. 

For the (VT) two-point functions, only J PC = 1 vector mesons contribute. In the LMD 
approximation, the leading short- distance behaviour then fixes everything but the mass of the 
lowest vector resonance, 

U L v r( P 2 ) = -(#)o^y. (4.7) 

It is quite remarkable that with this simple ansatz all the remaining leading short- distance 
constraints explicited in the previous section are met. In particular, for the quantities intro- 
duced in Eq. (|3.13|) we obtain 
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which satisfy ( 3.14 ). Note, however, that the LMD ansatze for T and Q in Eq. ( |4. 1|) do not 
correctly reproduce the subleading terms in the OPE for the correlator (VAP) given in Ref. |33 
The LMD ansatz for the invariant functions Tiy,A reads 

/LMD/ 2 2, , ^2^ _ W)o p 2 + q 2 + (p + q) 2 -c v _ N c M v 



2 ■ 



njLMUj 2 I ( I \2\ \yy/u f 1 ^ 1 \f 1 ^ J ~v 

n v {P,q,{P + q)) - ~ {p2 _ M , ){g2 _ M , ){p - q)2 , C v-^2 F 2 

njLMD / 2 _a fn , n ^ (^)p p 2 + q 2 -(p + q) 2 -c A N c M\ 



where cv,a are fixed by the Wess-Zumino-Witten anomaly [P] term. Again, these simple ansatze 
fulfill all the remaining leading short-distance requirements worked out in the preceding section. 
We note that for the vertex functions Tyy and r aa also the subleading terms in the OPE in 
Eq. ( |A.6| ) are reproduced by these ansatze. 

As already pointed out, it is sometimes necessary to generalize the ansatze for the invariant 
functions given above by including more than one resonance in each channel. This might be due 
to some additional constraints that are imposed on the Green's functions or in order to better 



reproduce experimental data involving resonances, as was argued in Ref. |32| . If we include, for 



instance, one additional vector resonance, the expressions for the invariant functions J-', Q, Hyr 



3 The LMD ansatz for the (VVP) three-point function was given in Refs. 
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and Hy read as follows (Ha remains, of course, unchanged) 

(V^)oP 2 [p 2 -q 2 - (p + qf] + Pf{p 2 , q 2 , (p + qf) 

2 ( P 2 -M2)(p2_ M 2)(g2_ M 2 )(p + g) 2 



^ + V,? 2 ,(P + <z) 2 ) = (Wo ^V + Pc( P 2 ,g 2 ,( P + g) 2 ) 



-.2 



^ MD+y (pw,(p+<7) 2 



where 



jtLM D+V / 2\ _ _ I-, - ,\ V +°VT 

UvT [P) ~ {W ° (V 2 ~ M 2 Vl )(p 2 ~ M 2 2 

(^)o P 2 q 2 [p 2 + q 2 + (p + q) 2 ] + P%(p 2 , q 2 , (p + q) 2 ) 



(p* - M 2 J(p^ - M 2 2 )(q* - M 2 ^ - M 2 2 )(p + qf ' 

(4.9) 



PF(p\q\(p + qf) = fip 2 + f2q 2 + f3(p + q) 2 + U, 
Pg(p 2 , q 2 , (p + <?) 2 ) = #iP 2 + 5 , 2? 2 + 5-3(p + g) 2 + #4, 
^(p 2 ,9 2 ,(p + g) 2 ) = /n(p 2 + g 2 ) 2 + /i 2 pV + /i 3 (p 2 + g 2 )(p + g) 2 + /i4(p + g) 4 

+/i 5 (p 2 + g 2 ) + h( P + q ) 2 + h 7 . 

The coefficients that appear in the polynomials Pp,Pc and P^ have to fulfill the following 
relations 

h + /s = -2cy T , 

h x + hz + h± = 2c VT , (4.10) 

in order to reproduce all short-distance constraints from the operator product expansion given 
in the previous section. Furthermore, the Wess-Zumino-Witten anomaly determines 

7 4tt 2 F 2 ■ 

Using the low-energy theorem that relates Tva to (VV — AA) one obtains, from the LMD+V 
ansatz for the latter correlator [45||, the relations 

A + / 2 = 2( gi + g 2 -M Vi -M 2 2 -M\ 

h = 2 [g, + M 2 Vi M 2 2 + M 2 A (M 2 Vi + M 2 2 ) - A ™f^ . ( 4 .11) 

The first of these conditions also guarantees that the next-to-leading short-distance behaviour 
of Tva is correctly reproduced. The vector form factor of the pion now reads 

2 2 

jpx.LMD+V, 2x _ , _g gig +94 19 x 
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Requiring that it behaves like 1/q 2 for large q 2 , leads to the relation 

9l = 2M 2 A . (4.13) 

On the other hand, in order to reproduce the subleading terms in the OPE for Typ, Eq. ( |A.6|) , 
with the LMD+V ansatz for Q yields a constraint which is independent from the previous one, 

g 1+ g 3 =-Ml. (4.14) 

Combining this equation with Eq. ( |4.13| ), we obtain the result 

2 „ 
g 3 = --M 2 A . 

Note that in contrast to the LMD case above, we can simultaneously fulfill the requirements 
from the asymptotic behaviour of the form factor Fy(q 2 ) and from the subleading terms in the 
OPE for the Typ vertex functions. 

Finally we note that the subleading terms in the OPE for the (VV|7r) vertex function Tyy in 
Eq. (|A.6|) are reproduced by the LMD+V ansatz for Hy without leading to further constraints 
on the coefficients hi. 

Let us briefly mention the ansatz for Tiy with one additional pseudoscalar resonance, dis- 
cussed in Ref. [ |32|1 , 

(Wo (p + qf [P 2 + q 2 + (p + qf] + P§ (p 2 , q 2 , (p + qf) 



H v MD+P {p\q\{p + q?) 



where 



2 (p 2 -M v )(q 2 -M v )((p + q) 2 -M 2 )(p + q) 2 ' 

(4.15) 



PH{p 2 ,q 2 ,{p + q) 2 ) = h^p 2 + q 2 ) + h 2 (p + qf + h 3 . 
The OPE leads to the condition 

h 1 = -M 2 , (4.16) 
and the Wess-Zumino-Witten anomaly yields 

t N c M v Mp 

h - 4^~Fjr- 

In contrast to the LMD case, the coefficients that appear in the invariant functions are 
no longer fixed unambiguously by the leading terms in the OPE. Here we have considered 
additional restrictions arising from the next-to-leading short- distance behaviour of the (V^A|7r), 
(VP\n) and (W^tt) vertex functions. Further information may be gained from the study of 
vertex functions like (VP\ai), (VV\p), etc. We shall not pursue this interesting line of thought 
here. Other sources of additional constraints might also be invoked, either from low-energy 
physics, from subleading terms in the OPE of the three-point functions or from processes 
involving resonances. We shall illustrate this point below. 

4 Actually, the addition of a single vector or pseudoscalar resonance is still not sufficient in order to reproduce 
the next-to-leading short-distance behaviour of (VAP) as given by Eq. (46) of Ref. J33[. 
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5 Determination of low-energy constants 



In this section, we shall use the ansatze of the previous section in order to obtain an evaluation 
of the low-energy constants involved in the chiral expansion of the three correlators under study. 
This is done upon performing the low-energy expansion of the invariant functions T LMD , Q LMD 
and Hy^ D and by subsequently matching it with the chiral expressions Q2.5| ) and (|2.6|). 

We start with the sector of even intrinsic parity, i.e. with the (VAP) correlator. For small 
momentum transfers, we can expand the resonance propagators 



P 



Ml 



Ml 



O 



and obtain (in these expressions, we have taken b = 2M 
after Eq. 



LAID 



(p 2 ,q 2 ,(p + q) 2 ) 

p2 



(p + qf 
(m\ 



M 2 



M V M\ \M 2 



M\ 



(5-13 



+ 



M 2 M 2 A 



Ml 



My, see the discussion 



2 Ml ) 2 M 2 M 2 + U ' 



■>LMD 



{p 2 ,q 2 ,{p + qf 



My 



My, 



My Ml, 



+ 0(p 4 ) 



q 2 (p + q) 2 

where 0(p 4 ) includes all possible higher order polynomials in p 2 , q 2 , (p + q) 2 . Comparison 
with the expressions ( |2.5| ) from ChPT yields the following solution (treating SU(2) and SU(3) 
together) 



and 



tLMD 
L 9 



tLMD 
^10 



-_Z. 



1LMD 
L 5 



LAID 



1 F 2 
1 ^0 



2 My 

1 F 2 



1 F 2 
1 r o 



4 My 



4 M 2 A ' 



r<LMD 
°78 



n LAlD 
°82 



riLMD 
°87 



„LMD 
-44 



LAID 



3F 2 



8 My 



+ 



F 2 



8 My Ml ' 



1 F 2 
1 r o 



F 2 



"17 



8 M£Mj ' 



c 



LAID 



1 F 2 



r,0 



8M£ 



+ 



F 2 



8 My Ml 



+ 



1 F 2 
1 r o 



8Mj ' 



LAID 



c, 



n LAID 
°89 



C 



1 F 2 
1 ^0 



51 



LAID 



4 My ' 

3F 2 



F 2 



-52 



4M£ 



2 .i/f .u-; ' 



LAID 
90 



C 



LAID 



53 



(5.2) 



(5.3) 
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The results for L 9 and L w agree with those obtained in Refs. after employing the two 

Weinberg sum rules [HDJ and using, in addition, the relation FyGy = Fq which follows from the 
assumption that the vector form factor of the pion Fy(q 2 ) satisfies an unsubtracted dispersion 
relation fEZU . 



Going through the same steps for the two correlators (VVP) and (AAP), we obtain the 
small momentum expansions 



n L v MD (p\ q \(p + q f) 



(p 2 ,q\(p + q) 2 ) 



(p + q) 2 
p 2 + q 



N, 



c 



8n 2 F 2 



2 ^~ 2 n Nr Ml 



My 

(Wo 

(p + g) 2 



8tt 2 F 2 



N, 



c 



24tt 2 F 2 



p 2 + q 2 ( 1 



iVc M 2 



24vr 2 F 2 



(g + g) 1 
M* 



(p + g)' 
Ml 



0(P 4 



+ C>(p 4 



from which we infer the following equations for some of the 0(p 6 ) low-energy constants A\ of 
Ref. |J in the odd intrinsic parity sector, see Eq. (|2.6|), 



.4 



LAID 



AA 



LAID 



j^LAID _|_ pLMD _ ^j^LML) 



o a LAID j_ A a LMD , a LAID a LAID 
0/± n -f / ±/l 17 "I" /i 2 3 — z/i 24 



x 21 



4^™ 



F 2 

8My 

8M V ' 
F 2 o 



N, 



c 



327T 2 My ' 



iV ( 



8M^ 96vr 2 Mi ' 
F 2 



\M% 



(5.4) 



The numerical values that follow from these expressions for the low-energy constants Cj and 



Ai are discussed in the next section. 



6 Comparison with the resonance Lagrangian approach 

In this section, we wish to compare the preceding determination of the low-energy constants 
with the approach which uses a Lagrangian with explicit resonance degrees of freedom. For 
definiteness, we use the resonance Lagrangian given in Ref. |Q which has often been employed 
in the literature to estimate the low-energy constants at order p 6 . In this Lagrangian a vector- 
field representation is used for the vector and axial-vector resonances. For convenience, we have 
written down this Lagrangian in Appendix |B[ Note that no pseudoscalar resonances appear in 
this Lagrangian. Furthermore, as stressed in Ref. p0 |, one has to add local terms from £ 4 , Q 



5 Where C n = X)fc+22=n ^(k.i) i n terms of the notation introduced in Eq. (1.1) 
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with fixed coefficients L\ es in order to correctly reproduce the QCD short- distance behaviour 
of certain Green's functions. We shall come back to this point below. 

The calculation of the (VAP), (VVP) and (AAP) three-point functions with the resonance 
Lagrangian yields again a result in agreement with the general solution of the Ward identities 
with the following expressions for the invariant functions (the coupling /3y is sometimes denoted 
by f x , see e.g. Ref. |§) 



F res (p\q 2 ,(p + qY) 



G res (p\q\(p + qY) 



and 



H r v ea (p 2 ,q\(p + qy 



nr(p\q 2 ,(p + qY 



p$(p + q y 



4L r g es + 4L[ e s 



(p 2 - M$) 



(q 2 - m\) 



p2q2 



(p 2 -M 2 )(g 2 
W>0 



f 2 - 2f v g v + 2V2 f v a v 
-f 2 A - 2V2 f A a A ) 

(-2f v fA(AW - A®)) 



M 2 A ) 



(6.1) 



F 2 q 2 (p + q) 2 
1 



AL r 9 es 



+ 



(p 2 - m 2 : 



(-2f v g v p 2 + 2V2f v a v q 2 - 4V2f v [3 v (p + q)' 



(q 2 



M\ 



2V2f A a A 
VvJaW 



(p 2 -M 2 )(g 2 



F$(p + q) 2 
p 2 



Ml) 



(A^q 2 -A^p 2 + 2B(p + qf 



(6.2) 



N, 



c 



(P 2 - M 2 \ 



(4V2f v h v ) + 



(q 2 - My] 



(AV2f v h v ) 



p 2 q 2 



(p 2 - M v )(q* _ M { 



v a v 



F$(p + q) 2 
p 2 



N, 



a 



24tt 2 



(P 2 



Ml] 



p 2 q 2 



(p* - M 2 A )(q2 - Ml 



(6.3) 



For momentum transfers that are small as compared to the resonance masses, we can ex- 
pand the propagators as sketched in Eq. ( |5.1|) . Therefore, the contributions from the resonance 
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Lagrangian start at 0(p 6 ) in the chiral expansion. The contributions A^ 2 \ A^ 3 \ B, and ay, <7a 
originate from the exchange of two resonances and start to contribute to the low-energy expan- 
sion at 0(p s ) only. 



6.1 Even intrinsic parity sector 

For the (VAP) correlator, we thus find the low-energy expansions 

(Wo 



F$(p + q) 2 
P 2 



M 2 



M\ 



(f v - 2f v9v + 2V2fi 



yOLy 



(-fl - 2v / 2/a«a) + 0(p 4 ) 



g res (p 2 ,q 2 ,(p + qf 



4L r 9 es 



F$q 2 (p + q) 2 

4- (-ZfvgvP 2 + 2V2f v a v q 2 - 4y/2 f v p v (p + qf 



Mi 

q 



M\ 



-2V2f A a A ) + 0{p 



Comparison with the expressions of the functions T and Q in ChPT, Eq. ( |2.5|) , leads to the 
following determination of the corresponding low-energy constants (again treating £77(2) and 
£77(3) together) 



/~tres 
°78 


_ res 
— t-44 


" AMy JV 


1 

+ 8 


M 2 ^ 


r~ires 
°82 


_ n res 
— C47 


1 1 


fv- 


- 1 1 / 


16 My 


16 My 


r~ires 
°87 


_ n res 
— c 50 


8My JV 


1 

~ 8 


M 2 A JAl 



2~^M fv(3v + 



AV2M 2 



fvPi 



T2~k fAaA 

1 



AV2M 2 a 



fACXA, 



Cres 
88 

Cres 
HQ 



rires 
°90 



1 1 



'51 



-52 



-53 



A My 



fv9v 



2My JV 

1 1 

'V^My 



11. 1 

rt/vSV H 



4M, 



72 M2 



V2M 2 



fv/3 



V 



(6.4) 



There are two main differences with the LMD ansatz of the previous section: the absence of a 
term (p + q) 2 in the low-momentum expansion of T res (j) 2 , q 2 ,{p + q) 2 ), whereas such a term is 
present in Q res (p 2 , q 2 ,{p + q) 2 ) but not in Q LMD (p 2 , q 2 ,{p + q) 2 )- 

In Table [l] (see e.g. Eq. (|5.3|) for the translation into the corresponding SU (2) constants q) 
we compare the numerical values for the low-energy constants in the LMD approximation 
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with those obtained from the resonance Lagrangian. As recalled in the introduction and as 
discussed in Ref. 



these numbers have to be understood as the values of the low-energy 
constants at the scale set by Ah, which we identify with the p mass My. We have used the 
values F = 92.4 MeV, M v = 769 MeV, and M A = 1230 MeV, as well as the two sets of input 
values for the resonance parameters listed in Table |2] as given in Ref. |33] (Set I), based on an 
ENJL model, and from Refs. [ 46|, [47H (Set II), extracted from resonance decays. If we allow 



for a relative error of about 30% in these values, a typical size for the uncertainty attached 
to a large-iVc estimate, the agreement is rather good, except for the cases of Css and Cgo- 
We postpone the discussion of some phenomenological implications of the differences shown by 
Table [I] to section 8 below. 



Table 1: Numerical values for the low-energy constants Cj (in units of 10 _4 /F 2 ) obtained from 
the LMD estimates in Eq. ( |5.3| ) and with the expressions derived from the resonance Lagrangian 
in Eq. (|6.4j ) for two different sets of input values for the resonance parameters. 





C*78 


C*82 


C*87 


C*88 


C*89 




LMD 


1.09 


-0.36 


0.40 


-0.52 


1.97 


0.0 


Set I 


1.09 


-0.29 


0.47 


-0.16 


2.29 


0.33 


Set II 


1.49 


-0.39 


0.65 


-0.14 


3.22 


0.51 



Table 2: Values for the parameters in the resonance Lagrangian from Ref. 
Refs. H, (Set II). 



(Set I) and from 





fv 


9v 


ay 


Pv = f x 


fA 


a A 


Set I 
Set II 


0.17 
0.20 


0.08 
0.09 


-0.015 
-0.014 


-0.019 
-0.025 


0.085 
0.10 


-0.0092 
-0.0067 



The values displayed in the second line of Table 1 were obtained by taking b = 1M\ in 
the LMD ansatz (ITT). If we had taken b = 4M|/3 instead, these values would have changed 



within an acceptable range: about 20% for C 78 and C 8 2, 16% in the case of C 8 g. The biggest 
variation occurs for C 8 s, around 30%, while C 8 7 and Cg remain unchanged, being insensitive 
to the value of b. Notice also that Lg is proportional to b, Lg = F^bj %M V M\. Changing b from 
2M\ to 4M|/3 decreases the value of Lg by as much as 33%, while leaving L w unaffected. This 
modifies the 0(p 4 ) prediction of the pion charge radius from (r 2 )y(b = 2M\) = 0.47 ± 0.13 
fm 2 to (r 2 ) v (b = 4M|/3) = 0.33 ± 0.09 fm 2 , somewhat lower than the experimental value 
(r 2 ) v = 0.439 ± 0.008 fm 2 ||]. 

At this stage, it is worthwhile to stress that it is not possible to find a one-to-one correspon- 
dence between the parameter sets of the resonance Lagrangian and of the LMD approximation 
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to large- Nq QCD. This directly reflects the differences mentioned above in the low-energy ex- 
pansions of the functions J^ LMD and Q LMD on the one hand, and of T res and Q res on the other 
hand. Nevertheless, some agreement can be obtained by adjusting the resonance parameters. 
When performing such a parametric comparison, one has to observe that the LMD ansatz al- 
ready encodes some additional information. For instance, the two Weinberg sum rules [ID] are 
fulfilled. In the LMD approximation they take the form F v = F 2 + F\ and F V M V = F\M\ 
and allow one to express Fy and Fa through the resonance masses and Fq. Furthermore, in 
the LMD approximation the identity FyGy = Fq holds. Finally, one has to make the following 
identifications between the parameters in the vector- and the tensor-field representation for the 



resonances: fy = Fy/My,gy = Gy/My and /a = Fa/ Ma, see Ref. |20 . 

Comparing the expressions for C90 and Css in the LMD approximation, Eq. ( |5.3| ), with 
those obtained from the resonance Lagrangian, Eq. (|6.4j) , and using FyGy = Fq we get 

Pv = 0. 

This removes the largest numerical discrepancies in Table |l|. Note in particular the huge 
cancellation in Cg| s for the values of (3y given in Table [2|. On the other hand, using the 
Weinberg sum rules and the identifications mentioned above one notices that 

/~iLMD — /~ires 
°87 — °87 • 

Thus we are left with the three equations C\ MD = C\ es ,% = 78,82,89, for the remaining two 
unknowns ay and a a- It turns out that this system of equations is inconsistent. We can solve, 
however, for ay, since in the difference CVs — Cs9 the term with /a a a drops out. From the 
requirement Cf 8 MD - C^ MD = - C 8 r 9 es we obtain 

T'2 



V2F M v (Ml + M v) 

Oly — 5 , — — U.U10 

8 Mi \\\ M y 



in remarkable agreement with the values quoted in Refs. 0, ^6], fllj , see Table ^|. On the other 
hand, requiring Cf 8 MD + AC^ 2 MD = Cf| s + 4Qf , leads to 



F « -0 



8 MyM\ 



This is not compatible with the spontaneous breakdown of chiral symmetry [I] and the Gold- 
stone theorem, which requires Fq 7^ 0. 

6.2 Odd intrinsic parity sector 

The low-energy expansion of the resonance expressions for the two correlators (VVP) and 
(AAP) gives the following estimates of the low-energy constants A{ 
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A r 2 es - 2A r 3 es - 4A r 4 es = 0, 
Ag + Ag-AAg = -V2f A h A 



Mi' 



3A{l s + 4A™ S + A r 2 f - 2A r 2 f - AA™ S = 0. (6.5) 

In Table |3| we compare the numerical values for the low-energy constants A* from the LMD 
estimates in Eq. (|5.4j ) with those from the resonance Lagrangian, Eq. ( |6.5| ). We have introduced 
the following notations for the combinations of low-energy constants that appear in (VVP) and 
(AAP) 

A VtP 2 = A 2 -AA 3} 
Av,(p+q) 2 — A 2 - 2A 3 - 4A 4 , 

A AtP 2 = A u + A 23 - AA 2i , 
A Mp+q) 2 = 3A U + AA 17 + A 23 - 2A 24 - AA 25 . (6.6) 



Table 3: Numerical values for the combinations of low-energy constants A t defined in Eq. ( |6.6| ), 
in units of 10 _4 /F 2 , obtained from the LMD estimates in Eq. (|5.4| ) and from the resonance 
Lagrangian in Eq. ( |6.5D (Set I). 





Ay p 2 




Aa, p 2 


A A ,(p+q) 2 


LMD 
Set I 


-1.11 
-1.13 


-0.26 
0.0 


-0.14 
-0.096 


0.040 
0.0 



The agreement for the low-energy constants Ay t A,p 2 is quite good, whereas the two ap- 
proaches give different results for Ay t A,(p+q) 2 - in the expressions for the low-energy constants 
Ai from the resonance Lagrangian ( |6.5|) we used the ENJL estimates (Set I) 



fv h v = ^T-2^ 1 + 9a) = 0.0055, 
16^^ 8 

f A h A = -^~-g A {l + g A )= 0.0012, (6.7) 
with Nq = 3 and g A = 0.65, see Ref. p8[, in particular the Erratum. We shall discuss estimates 



for the low-energy constants Ay )P 2 and Ay^ p+q y beyond the LMD approximation in section [873 . 

Again, it is impossible to find a one-to-one relation between the parameters of the resonance 
Lagrangian and those that describe the LMD ansatz. The reason lies in the absence of a 
term proportional to (p + q) 2 in the low-energy expansions of both Tiy s and 7i r A es . We notice 
nevertheless that requiring A V K /P = Ay^ p2 and Aj&P = A r A e ^ p2 leads to the relations 

1 / AT F 2 \ 

f v h v = — —°- = 0.0054 , 

JV V y/2 V 327T 2 8M 2 , J 
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fAhj 



l_ (Ni 
V2 \96 



c 



8M 2 A 



0.0017. 



in rather good agreement with Eq. ( |6.7| ). 

We conclude that, although an adequate adjustment of the parameters of the resonance 
Lagrangian of Ref. [[38] can bring the determinations of the low-energy constants considered 
here to a reasonable numerical agreement, there is no way to establish an algebraic equivalence 
between the two approaches. We have focused here on the Lagrangian of Ref. which is the 
most complete as far as couplings among resonances and to the Goldstone bosons are concerned. 
Therefore similar conclusions will also hold for other existing Lagrangians with resonance fields, 
see Refs. 



14, 15, 16, 17, 18, 19, 20] and references therein. 



7 QCD short-distance constraints on the resonance La- 
grangian 

The discrepancies between the estimates for the low-energy constants at 0(p 6 ) from the LMD 
ansatz and from the resonance Lagrangian, as observed in the previous section, can be traced 
back to the different high-energy behaviours of the corresponding Green's functions in the two 
approaches. In fact, as we shall show in this section, the Green's function derived from the 
resonance Lagrangian are incompatible with the QCD short-distance constraints. 



7.1 Two large momenta 



We first consider the limit when the two momenta become simultaneously large, see Eqs. (3.2) 
- Q3.4j) . The invariant function T res from Eq. ( |6.1[) behaves as 

hm F es ((\p)\ (\q) 2 , (Xp + \q) 2 ) = .J^, , 2 Ul 9 + 4L 10 + r FV + r FA + r FVA 



A 4 Fn 



(r F y + r FVA )M v (r FA + r FVA )M\ 
p 2 (p + q) 2 q 2 (p + q) 2 



where r FV , r FA , and r FVA denote the coefficients of the terms p 2 / (p 2 — My), q 2 / (q 2 — M A ), and 
(p 2 q 2 )/{(p 2 — M v )(q 2 — M A )) in Eq. Q6.1]) , respectively. Agreement with the OPE result from 
Eq. ( |3. 2D at order 1/ A 2 can be achieved if the constraint 

4L 9 + 4L 10 + f 2 - 2f v9v + 2V2f v a v - f \ - 2v / 2/a«a - 2f v f A (A^ - A^) = (7.2) 

is fulfilled. However, at order 1/A 4 , we observe that T res is not compatible with the OPE result 
in Eq. ( |3.2| ), since the term ~ ^-/{p 2 q 2 ) is missing in Eq. ( |7. 1| ) . 

As was shown in Ref. p0| , requiring agreement of certain Green's functions with the short- 
distance properties of QCD uniquely determines the low-energy constants Li at OQo 4 ), even 
though the contributions of the resonance Lagrangian in the vector-field representation only 



21 



start at (9(p 6 ). 6 In fact, one has to add local terms from £4 to the resonance Lagrangian in 
order to obtain the correct short-distance behaviour of these Green's functions. Analogously, 
at 0(p 6 ), one might try to add local counterterms from Cq. This is, however, not enough to 
bring the function T res in agreement with the OPE. One has as well to add new terms involving 
both resonance fields and additional derivatives. A similar observation, concerning the (VVS) 
correlator, was made in Ref. [31]: the short-distance behaviour of the three-point function 



{VVS) becomes consistent with the OPE only if a term ((V^ — (/y/2) /+) V S) is added to 
the resonance Lagrangian. In the present case the situation is more involved, i.e. several new 
terms would have to be added. We have not undertaken the task to construct them explicitly. 
One can show, however, that if one matches T res with the constraints imposed by the OPE, the 
local counterterms at 0(p 6 ) have to be adjusted in such a way that one finally obtains the same 
values for the low-energy constants as with the LMD approach. We caution the reader that the 
mere fact of using a tensor field representation for the resonances does not, by itself, guarantee 
to yield the same estimates for the low-energy constants as with the LMD ansatz. Actually, it 
was pointed out in Ref. |33| that the resonance Lagrangian with a tensor field representation 
leads to a correlator (VAP) (49| which does not have the correct short-distance properties. 

We may perform a similar analysis for Q res from Eq. ( |6.2j ). In this case, one may obtain 
agreement with the QCD result of Eq. ( |3.3| ), provided the coupling constants in the resonance 
Lagrangian are adjusted as follows 

f v g v M 2 + V2f v a v M 2 A - V2f A a A M v = ^ , L 9 = ^f v9v , 

A^ = y/2^, A^ = V2^ } (3 V = 0, 5 = 0. (7.3) 

JA JV 

We note, however, that in deriving these constraints we have not taken into account possible 
contributions to Q res from the new local terms that have to be added to the resonance La- 
grangian in order to make T res compatible with the OPE. The expression for L 9 in Eq. (|7.3| ) 
agrees with Refs. ]T8| , E0|, whereas the relations for A^ and A^ agree with those given in 
Ref. 



The short-distance behaviour of the functions TCy s and Ti r A es in Eq. ( |6.3|) in the odd intrinsic 
parity sector is given by 



lmr^((A.) 2 ,(A#,(A, + A#) = " A8f?(p + ff)a 



V +8V2f v h v -4fta v 



57T 



(#>oM, 



A 4 F 2 



p 2 (p + q) 2 q 2 (p + q) 2 



'A 6 " 

6 If a tensor-field formulation is used for the vector and axial-vector resonances, the contributions from the 
resonances start already at order p and lead directly to the usual estimates for the Li. On the other hand, 
several couplings which appear in the resonance Lagrangian of Ref. [ ]38| , like ay , fiv or a a , cannot be written 
down in the tensor-field representation, at least not without introducing additional derivatives. 
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lim nr((\p) 2 ,(M) 2 ,(\p + MY 

A— >oo 



\ 2 _F 2 ( P + q) 2 

(#}om; 



N c 

247T 2 



+ 8V2 f A h A - Af 2 A a A 



\*F 2 



4V2f A h A -4fla A | AV2f A h A -Af 2 A a A 
p 2 (p + q) 2 q 2 (p + q) 2 



+0(±) . (7.4) 

Comparison of these expressions with the OPE result of Eq. ( |3.4| ) leads, at order 1/A 2 , to the 
constraints 



+ 8V2f v h v -Af 2 a v = 0, 



8tt 2 

j + 8V2 f A h A -4f 2 A a A = 0. (7.5) 
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However, at order 1/A 4 , we again observe that 7iy s and 7i r A es are not consistent with the 
OPE result, since the terms ~ ^/{p 2 q 2 ) are missing in Eq. ( |7.4| ). 

7.2 One large momentum 

In section ||] we have also derived the short-distance properties of the Green's functions 
(VAP), (VVP) and (AAP) when the relative distance between only two of the currents becomes 
small. In certain physical applications only this limit is relevant and we shall now discuss the 
corresponding constraints on the parameters in the resonance Lagrangian. We note, however, 
that it is not possible to satisfy simultaneously all the constraints given below. Moreover, 
some of these constraints are in contradiction with the relations derived in the previous section. 
These inconsistencies can again be traced back to the fact that the Green's functions derived 
from the resonance Lagrangian do not correctly reproduce the QCD short-distance behaviour. 

The constraints on the functions T and Q from Eqs. ( |3.13| ) and fl3.14| ), when the space- 
time arguments of the vector and axial-vector currents coincide, can be satisfied provided the 
resonance parameters obey the relations 

A^ = V2^, B = -V2^. (7.6) 

J A J A 

Furthermore, we recover the usual resonance estimate L10 = — /y/4 + /j|/4 and the first Wein- 
berg sum rule M 2 fy = Fq + M A f A . Imposing the relation L 9 = fygy/2, we furthermore get 
the relations 

A m = ^2^, M 2 f 2 - M 2 f v g v + V2M 2 A f v a v + V2M 2 f A a A = . (7.7) 
Jv 

The second limit, when the distance between the vector current and the pseudoscalar density 
becomes small, corresponding to the constraints ( |3.15| ), can be satisfied provided 

A (3) = v /2^i, /3 y = 0, B = 0. (7.8) 
Jv 
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We also get the result L 9 = f v g v /2. Imposing the usual resonance estimate for L 10 , we obtain 
the additional relations 



yCty = — : 



J A 



(7.9) 



Finally, the constraints (|3.16| ), when the arguments of the axial- vector current and the 



pseudoscalar density coincide, lead again to the relations (|7.6|) . Using in addition the usual 
resonance estimates for Lg and Lio, we obtain 



A^ = V2^, fy 9 y-^-V2f_ 



A&A 



Y 



F 2 
2M V 



(7.10) 



We note that we recover in the first case, corresponding to Eqs. ( |3.13|) and (|3.14| ), all the 
constraints from the leading terms in the expansion in 1/A, when all momenta in T and Q 
become large. In the latter two cases we get, however, only a subset of these relations. 

The equations ( [7.6|) - (|7.10| ) arising from the three different limits are perfectly compatible. 
If however we combine them with the first relation in Eq. ( |7.3[ ) from the OPE constraint when 
all momenta in Q become large, we find a contradiction. 

We now turn to the functions (VVP) and (AAP). We first consider the case when the 
two vector currents are taken at the same point, which is relevant, for instance, in the decay 
P — > where P = tt°,t] (see the discussion in Ref. [43|) or when the space-time arguments 



of the two axial-vector currents in (AAP) coincide. The corresponding constraint ( |3.17| ) can 
be satisfied, provided the resonance parameters fulfill the relations 



V2f x h 



x 



c 



1 F 2 



32tt j 



8M 2 ' 



fx^x = c x 



N, 



c 



1 F 2 



32vr 2 



4Mi' 



X = V,A, 



(7.11) 



with cy — 1, ca — 1/3. In the case when the space-time argument of one of the vector currents 
and of the pseudoscalar density in (VVP) coincide, the constraint from Eq. ( |3.18| ) can be 
satisfied, if 

N c 1 F 2 2 N c 



V2f v h 



V 



III 

A My' 



1 F 2 

2M 2} 



32vr 2 A My 1 JV 32tt 2 

where we have also used the short- distance constraint ( |3.5|) on Uvt(p 2 
(AAP), see Eq. ( ^.19| ), yields the relations 



V2f A h A 



c 



967T 2 ' 



N, 



o 



96tt 5 



(7.12) 
The same limit in 

(7.13) 



We recover in all cases the constraints ( [7.5| ) when all momenta in 7iy : A become large. We note, 
however, that although one can individually satisfy the constraints from the OPE, the relations 
( |7.12| ) and ( [7.1 3| ) are incompatible with Eq. ( |7.11| ). 
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8 Resonance contributions to pion form factors 



In this section, we discuss a few phenomenological applications of the various ansatze consid- 
ered previously. The first two examples, the vector form factor of the pion and the radiative 
pion decay, involve the low-energy constants C$ that were determined in section 5. The third 
example, the pion-photon-photon transition form factor, illustrates a situation where the MHA 
goes beyond the simplest LMD approximation. 



8.1 Vector form factor of the pion 

There are two combinations of renormalized low-energy constants c[ from the chiral Lagrangian 



Cq that enter in the vector form factor of the pion, see Ref. J50] 



r vl = -16c£ - 4^ 5 - 8c r 53 , 
r r V 2 = -44! + 4c r 53 . 

The resonance estimates for r vl and r V2 given in Ref. [JST]] read (recall that f x = fly) 

2V2fJ V F 2 



r r vi S (.H = M v) 



gvfvF 2 



r v T^ = M v ) = (8.1) 



Using the resonance estimates for the constants c[ es given in Eqs. ( |6.4| ) we obtain the same 
result for ry r 2 ei ', if we identify the pion decay constant F with F . 
With the LMD estimates from Eqs. ( |5.3| ) we get 

H 

Mi 



r v L 2 MD (» = M v ) = ^, (8.2) 



which agrees with the result given in Eq. ( ^.1| ), if we use the relation FyGy = Fq which is 
valid within the LMD ansatz. From the present analysis, we cannot obtain an estimate for r vi , 
which describes the quark mass corrections to the value of the vector form factor at vanishing 



momentum transfer, see Ref. 51 



8.2 The decay it — ► ez/ e 7 

The decay n(p) — > ev e ^(q) is described by two form factors V and A. The contribution from 
£ 6 to A can be written as [5(| 

A((p- q f) = My A1 + ( P - q y A2 , 



with 



Al 



' A2 



48c r e - 16c r M + 8c r 35 - 8cL + 16c^ - 16c r i7 + 8c 



-44 



-16 



-47 



50 ' 



8c 



li 



16c r 50 



+ 4c: 



51 
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With the estimates for the low-energy constants c\ obtained from the resonance Lagrangian, 
Eqs. (|6.4|), we obtain 



A2 



(jj, = M v ) = F 2 j^ (f A + f A a A 2y/2) « 0.55 • 10" 



13) 



where we used the same values f A = 0.080, a:^ = —6.66 • 10 , as in Ref. [52] and F ~ F ( 



92.4 MeV, and M A = 1230 MeV. 

On the other hand, using the LMD estimates given in Eqs ( |5.3| ) we obtain 



r A2 (t* = Ml 



V 



M V M\ 



F 4 
2^° 



0.18 • 10" 



5.4) 



The LMD estimate differs by a factor three from the value given in Eq. (|8.3|) . In fact, since 
we obtain a different relative sign in our result for r^ 6 * as compared to the expression given 
in Ref. [521, the LMD estimate is even a factor of five smaller than the value T A ^ es (^i = My) ~ 
0.89 • 10~ 4 used in that paper. 

Since the same combination of resonance parameters that determines r r A2 also enters in the 
decay amplitude for a\ — > 7c-f, the LMD ansatz predicts a decay rate more than an order of 
magnitude smaller than the usual value |k|. However, the experimental situation concerning 
this decay is far from being clear f54|, see the remarks following Eq. (62) in Ref. [|33|. 



8.3 Pion-photon-photon transition form factor 

Up to now, we have shown that a very minimal ansatz allows to take into account the leading 
asymptotic behaviours of the three point functions (VAP), (VVP) and (AAP). This simple 
representation would in general certainly not be sufficient if additional information were added. 
However, within the large-iV^ framework considered here, one always has the freedom to go 
beyond the LMD representation and to add other zero-width resonance states. In this section, 
we wish to illustrate this point by considering the form factor •F 7r7 * 7 *((7i, (?f) that describes the 
transition between a pion and two (possibly off-shell) photons in the chiral limit. This form 
factor is defined as 

• J d 4 xe 1 ^ (OlTO^^O)}^)) = e Fa //^ ff ( 9 2 , (p - qf) (8.5) 

or still 



For both photons on-shell, the value of this form factor is fixed by the Wess-Zumino-Witten 
anomaly term, 

N c 

^,(0,0) = -^-. 
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Many studies have been devoted to this form factor in the past, see j|2[ [55| and references 
therein. In particular, its behaviour in the limit Q 2 — > oo, uj fixed, with — Q 2 = (q\ + gf)> 
uj = (qf — gf)/(g 2 + si); has been investigated, with the result 



2 2 4F 



(8.6) 



3 g 2 

for — 1 < u; < 1. The function f(u>) can be expressed in terms of the pion distribution function 



du 



(1-u)(1 + uj) +u(l -u) 



[1 + 0{a s 



normalized as 



du (p n (u) 



1 . 



This last condition is sufficient in order to study the limit — q\ = —q\ 
to the result [KB 



Q 2 /2 — ► oo and leads 



J~ 777*7* ( „ , 



Q 2 

2 



g 2 

2 



4F n 1 



1 



3 Q 



-^ + o — 



Q 1 



The function ip n (u) is only known asymptotically, and this asymptotic expression is reliable 
only for u not too large [57|, e.g. |u;| < 1/2. The case of one on-shell photon corresponds to 
\u\ = 1, so that the coefficient of the l/Q 2 fall-off of the form factor T^*^» (— Q 2 , 0) is actually 
not known. Depending on the assumptions made or the ansatze considered for cp^lu), different 
results have been obtained in the literature. 

The LMD ansatz for Tiy^q 2 , gf, (g x + g2) 2 ) reproduces these results for \u\ < 1 with 



./ 



LMD, 



UJ 



1/(1 — uj 2 ). On the other hand, taking q\ = and letting Q 2 = —q\ become 



large, we obtain 



LMD i 

7T7*7* \ 



-Q 2 , 0) ~ const. 



In order to recover the l/Q 2 behaviour for 



\0J\ 



1, we need to go beyond the LMD approx- 



imation and add a second vector resonance (adding a pseudoscalar resonance f32| would not 
help to improve the situation in the present case, since the l/Q 2 behaviour forces hi = in the 
ansatz for Tiy, in contradiction with the relation ( 4.16[) ). From Eq. ( f4.9|) we obtain 



LMD+V 



F q\q\ \q\ + ql] + h x {q\ + q\f + h 2 q 2 q 2 + h 5 (q 2 + q 2 ) + h 7 



(qI-M 2 Vi )(q 2 -M 2 2 )(q 2 -M 2 Vi )(q : 



5.7) 



The behaviour of F^® +v (q 2 , q 2 ) for Q 2 large and uj fixed, with |u;| < 1, is the same as in the 
case of T^l®{q\, g|). However, if we now set q\ = 0, we obtain 



-pLMD+V, 

7T7*7* ' 



-g 2 ,o) 



hQ 4 - h 5 Q 2 + h 



3 M v M 2 2 (Q 2 + M Vi )(Q 2 + M 2 2 : 
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Imposing that this expression exhibits the 1/Q 2 behaviour for large Q 2 requires that hi vanishes, 
which then gives 



In the absence of a reliable prediction for the coefficient that governs the 1/Q 2 behaviour, this 
still leaves the parameter h 5 undetermined. Additional information may be obtained from the 
fact that the form factor ^v r *- r *(g r i, ?|) can a l so be related to the decay p + — > 7r + 7, whose 
amplitude is given by 



( qj - M? 

FyMy 



A(p + ^n + 1 )= - lira hm ™ S'S* ' F^Ati, <£) ■ (8.10) 



Note that we obtain the same relation for A(p° — > 7r°7). In the latter case, however, p — uj 
mixing would have to be taken into account for a realistic calculation. Furthermore, we have 
defined the coupling Fy of the p meson to the vector current by 

{0\V^)\p\p)) = 5 ah FyM v e^ 

where denotes the polarization vector of the p meson. From Eq. ( |S.10| ) we obtain (see also 
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(ItFy \ A{p + -> 7T+ 7 ) \im q 2 M 2 2 v(ql - M 2 )T^ rT (qj, qf) 

1 1 - ' ' — 1 + X . (o.ll 



\ M v J A(tt° -> 77) hm ? 2^o, 3 |^o(gi ~ M^)^ rr (qf,q^ 



The observed value T = 68 ± 7 keV for the decay width p + 7r + 7 |53| yields x = 0.022 ± 0.051 
[0. The LMD ansatz for leads to x LMD = -(4vr 2 /N C )(F 2 /M v ) = -0.19, far from the 

experimental value f58|. Starting instead from the form factor jF 7 ^f^ +y leads to 



1 + x) LM D+v - t\ IfTJTJT^ I 1 



4tt 2 F 







1-M 2 JM 2 2 )\ N c M Vi 



M 2 x h x h 5 



M 2 2 M 2 2 M V2 



Setting hi = and solving for h$ gives h 5 = 6.3 ± 0.9 GeV 4 , where we have taken My 1 = 
769 MeV, My 2 = 1465 MeV for the resonance masses and F = 92 A MeV. 

Actually, the form factor ^ r 7r7 * 7 *(— Q 2 , 0) has been measured in the space-like region by the 
CLEO collaboration, see Table 1 in Ref. p§, over a wide range of Q 2 , 1.5 GeV 2 < Q 2 < 9 GeV 2 . 



A fit of the expression T^l^ v {— Q 2 , 0) with hi — to these data yields 

h 5 = 6.93 ± 0.26 GeV 4 , (8.12) 

with x 2 /dof = 7.00/14 = 0.50. Keeping also hi as a free parameter yields instead 

hi = -0.01 ±0.16 GeV 2 , 
h s = 6.88 ± 0.61 GeV 4 , 
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with x 2 /dof = 6.99/13 = 0.54. The results for /i 5 from both fits are compatible with the value 
extracted above from the decay p + — > 7r + 7. 

Finally, one might ask to which extent the inclusion of a second vector resonance into the 
ansatz for Tiy modifies the determination of the corresponding combinations of low-energy 
constants. The only combination which can be fixed from the knowledge of J r 7r7 * 7 * {—Q 2 , 0) 
alone is Ay >p 2 (see Eq. ( |6.6|) ), 



a lmd + v _ HtJ^ _ N^J_ ( 

A ^ 2 - 8 M* M V2 32vr 2 M^y + M*J> [ * A6) 



since the other combination, A V( J^\ 2 , involves h e . With the value of h 5 obtained in Eq. ( |3.12|) , 



we find, in units of 10~ 4 /F^ , A V ^ 2 D+V = -1.36, i.e. about 20 % away from our LMD estimate 
reported in Table |3|. The difference is well within the 30% relative error that we attribute to 
the approximations considered there. 

Another approach was followed in Ref. [32], where an additional pseudoscalar resonance it' 



was included in the ansatz for Tiy that satisfies the OPE constraints (LMD+P, see Eq. (|4.15| )). 
As noted above, this ansatz will, however, not correctly reproduce the 1/Q 2 behaviour of 
< F Brr (_Q2 j o) at large Q 2 . In this case, A v<p 2 = -N c {l+x)/32ix 2 M v and from Eq. (|8~TT|) one 
obtains the result, again in units of 10- 4 /F 2 , A V ™ 2 D+P = -1.40, close to our LMD+V estimate 
from the fit to the CLEO data. In a similar way, the low-energy constant Ay^a receives 
within the LMD+P ansatz an additional contribution proportional to the decay amplitude 



A(n' — ► 77) [52]]. Since this decay has not yet been observed experimentally, our LMD estimate 
for Av,(p+q)n from Table ^| is probably not strongly modified by the addition of a pseudoscalar 
resonance. 

The low-energy constant Ay tP 2 also describes the contributions from the counterterms at 
order p 6 to the slope b n of the form factor at the origin 



A(tt° 77*(g 2 )) dq 




At 0(p 6 ) this slope also receives contributions from chiral loops, therefore, b n = b l ° ops + b^ T , 
with b^ T = —327i 2 Ay jP 2/N c . The Particle Data Group gives the value a n = M 2 b m = 
0.032 ± 0.004 ]S§. From our LMD+V estimate above we obtain b% T > LMD+v = 1.67 GeV" 2 
or, equivalently, a c T > LMD + v = 0.031. To this value, one should add the contribution from the 
chiral logarithms, evaluated at the scale \i ~ M p , a l ° ops ~ 0.005 |60[| , which represents a 20% 



effect. We note that the value for a n used by the PDG is essentially the one reported by the 



CELLO collaboration |61j| . In the latter paper a simple VMD-inspired pole ansatz was fitted 
with their data for the form factor in the space-like region for 0.5 GeV 2 < Q 2 < 2.7 GeV 2 , not 
taking into account contributions from Goldstone boson loops at low Q 2 . 
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9 Conclusions 



In this article we have studied the QCD three-point functions (VAP), (VVP) and (AAP) in 
the three-flavour chiral limit in order to obtain resonance estimates for some of the low-energy 
constants that appear at 0(p 6 ) in chiral perturbation theory in the meson sector (even and odd 
intrinsic parity). We have compared the results that have been obtained in the literature using 



a Lagrangian that includes resonance fields |38] with those evaluated within the framework of 



an approximation of large-iVc QCD combined with information on short-distance properties. 
In certain cases, we have found substantially different results for the estimates of the low- 
energy constants obtained with these two methods. We have pointed out that this is due to 
the fact that the Green's functions derived from the resonance Lagrangian do not correctly 
reproduce the QCD short-distance behaviour. This defect can be repaired, but at the expense 
of introducing, into the resonance Lagrangian, certain local contributions. The difference with 
the similar situation at the 0(p A ) level lies in the fact that these local contributions cannot be 
restricted to terms already present in the 0(p 6 ) chiral Lagrangian, but also involve terms with 
resonance fields and higher order derivatives. This feature, already noticed in a particular case 
in | 3"I|| , seems to be of a general character. A general construction remains to be done, and 
appears to be a much more complicated task than at 0(p A ). 

We note that although in general the short-distance behaviour of the Green's functions 



derived from the resonance Lagrangian [38] is incompatible with QCD, one can sometimes 
reproduce the results from the operator product expansion, by adjusting the resonance param- 
eters accordingly. There are, however, certain cases where this is not possible. The numerical 
values for the low-energy constants can then be very different. In particular, whereas both 
methods lead to identical estimates for the resonance contributions at order p 6 in the vector 
form factor of the pion Fy(q 2 ), our estimate for the resonance contribution in one of the form 
factors for the decay tt — ► ez/ e 7 is a factor of five smaller than the results quoted in Ref. p>2] . 



Of course one might argue that the short-distance behaviour of Green's functions, i.e. their 
behaviour at very high energies, is irrelevant for the determination of low-energy constants in 
chiral Lagrangians starting from a resonance Lagrangian that is supposed to be valid only in 
the intermediate energy region anyway. We think, however, that taking into account the QCD 
short-distance constraints is a good guiding principle to avoid (some of) the ambiguities when 
working with resonance fields, as was shown at order p A in Ref. ||20|| . Moreover, in certain cases 
one needs also integrals of these Green's functions, for instance, to estimate the low-energy 



constants that appear if virtual photons |34J] or leptons [|37| are included in chiral perturbation 
theory, see the discussion in Refs. [B31, The case of the counterterms of the effective 



Lagrangians |)2| in the |A5| = 1 or \AS\ = 2 sectors of the Standard Model presents very 
similar features p3] , |5^[ . In these applications it is crucial that the Green's functions respect 
the short- distance constraints in order to obtain ultraviolet finite results or to implement a 
correct matching between long and short distances. 

Finally, we wish to add a few remarks concerning the methodology followed in the present 
study. We shall not come back on the use of the large-iVc framework, which seems to be 
unavoidable once correlators of higher rank than two-point functions need to be considered. 
Our approach deviates however from a full large-iVc limit of QCD by at least two aspects. 
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For one thing, we have only imposed constraints coming from the leading (for the three-point 
functions) and next-to-leading (for the vertex functions) short-distance properties of QCD. We 
have considered neither the effects of higher dimension operators in the Wilson expansions, nor 
have we included QCD corrections to the short-distance behaviour (recall however, that for 
the cases treated here, the corresponding Wilson coefficients had no anomalous dimensions). 
For the other thing, we have truncated the mesonic spectrum of large-iV<7 QCD to the mini- 
mal number of resonances necessary in order to fulfill the short-distance constraints that were 
considered. Both approximations are to a large extent interdependent. It is clear that, say, 
the simplest LMD ansatz will at some point fail to reproduce the subleading short- distance 
behaviour. On the other hand, the knowledge of the subdominant operators in the OPE, or 
other constraints, will fix additional parameters that have to be introduced if one goes beyond 
the LMD approximation. In this sense, the framework within which we have been working is, 
given the necessary amount of work, improvable. As an illustration, we have considered the 
pion-photon-photon form factor F„- 7 * 7 *(— Q 2 , 0). The experimentally observed |Kj l/Q 2 fall-off 
of the form factor cannot be reproduced with our LMD ansatz for the invariant function TCy. 
Including one additional vector-resonance, we obtain a representation for this form factor that 
fulfills all constraints from the leading terms in the operator product expansion and that fits 
the experimental results successfully if we adjust some of the unknown parameters that enter in 
the generalized ansatz for TCy- We note that this is not possible if one includes a pseudoscalar 
resonance instead. This new ansatz, with the phenomenologically determined parameters, is 
furthermore compatible with the observed decay rate for p + — > 7r + 7. Finally, the estimate for 
one combination of low-energy constants Ai in the odd-intrinsic parity sector changes by about 
20 % with this new ansatz. This difference is within the 30 % relative error that we attribute 
to the approximations considered here. The determination of the low-energy constants seems 
thus to be stable against the inclusion of higher mass resonances. 

It is clear that the three-point functions we have considered do not, by far, exhaust the 
whole set of low-energy constants of the chiral Lagrangian at the 0(p 6 ) level. Other three point 
functions, as well as higher correlators, need to be studied in a similar way for that purpose. 
Some of them (in particular, those describing the low-energy constants related to quark mass 
corrections) will include the scalar densities. Recent studies have emphasized that in the ++ 
channel, predictions relying on the large- Nc picture might not be very reliable because 
of the strong 7T7T interaction in the S-wave. It has been suggested that a more appropriate 
treatment would rather require to consider the limit where both Nc and Np, the number of 
light flavours, become large, with a fixed ratio Np/Nc fl66| . Also, in the case of two- and three- 
point functions, the quantum numbers of the resonances that can contribute in the large- de- 
limit are entirely fixed by the quantum numbers of the quark bilinears involved. This is no 
longer true for n-point functions with n > 4 [21|. We leave these and other interesting issues 
for future work. 
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A Operator product expansion for vertex functions 

In this appendix we discuss the short-distance behaviour of the vertex functions 



( r vA)% c (q,p) 



J d*xe^(0\T{V«(x)A b u (0)}\n c (p)} , 

: /^(o|t{i/;(x)p 6 (o)}|/( p )), 

= J d*xe^(0\T{V«(x)V u b (0)}\7r%p)) , 
(J uCfop) = fd 4 xe^{0\T{A;(x)Al(0)}\n c (p)) . 



( r w)t;(q,p) 



(A.l) 



They are related, according to the LSZ reduction formula, to the (VAP), (VVP) and (AAP) 
three-point functions as follows: 



( r VA)lt c {q,p) 



rabc 



W>0 



r 



/T/X f (gp ~ q)»(p - q)u 

m ° { — ^ 



F(q 2 ,q-p) 



+PtxM, v - q)F(q 2 , q-p) + Q»Aq, v - q)G(q 2 , q ■ p) 

lim p 2 F(q 2 ,(p-q) 2 ,p 2 ) , 

p^^O 



G(q,q-p) = limp Q{q ,(p-q) ,p ), 

p 2 ^0 



(A.2) 



and 



(XvpT^ (q,p) 



rabc 



'~' JA { ff^ + {q2p " " {q ' p)q ^ iq2 > q ■ p) 



■W)o- 



Q(q,q-p) = \imp G(q ,p ,(q-p) ), 

p^^O 



(A.3) 



since the invariant function ^ r (q ,2 ,p 2 , (<? — p) 2 ) does not contain a 1/p 2 Goldstone boson pole. 
Finally, 



( r vv) a ^(q,p) 



t J^e» l/aP q a pPd abc H v (q 2 ,q-p), 
(Wo 
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K x (q ,q-p) 



i^^e^ q a p^d ahc n A (q 2 ,q-p), 
(W)o 

lim p 2 H x (q 2 ,(p-q) 2 ,p 2 ), X = V, A . 
p — >o 



(A.4) 



with 7ix(q 2 , q • V) — Hx{.q 2 — 2g • p, — q ■ p), as a consequence of Eq. ( ^4| 
For general Dirac matrices r 12 we obtain the OPE 

lim / dW^)-(0|T{(^r 1 ^V')^)(^r 2 ^)(0)}|7r c (p)) 



^/^^(ro/^J (r 



'A^ x 



XL 



X 
+ 



iS^y^Ol : V?(x)^,l(0) : |tt c (p)> + iS(-x) 5 a S IL (0\ : t/^W/vO) : \n c (p)) 



(T S 



/J 



XL 



9 



5^(01 : ^(0)^(0) : |7r c (p)> 



5 IL (0\: ^(0)V/3,j(0) :|tt c ( P )) 



<9 



A fc 



7 



) 7 I 2 



XL 



^(0| : (D^)?(0)^,l(0) : |vr c (p)> 



- 4 «/i(0hfi(O)(W):|/(p)) 



1 



(A.5) 



up to possible 0(a s ) corrections. Note that the indices I, J, . . . label both flavour and colour 
and that in this Appendix (X a /2)u denote the Gell-Mann matrices in flavour space and the 
unit matrix in colour space, with tr(A a A b ) = 2NcS ab . From invariance under parity, Lorentz, 
flavour SU(3)v and colour SU(3)c transformations the matrix elements involved in Eq. (|A.5 ) 
can be expressed as 



: ^(0)^r(0) : |7r c (p)> 



E 



LI 



(il 5 ) s a (0\ : $i75^-VO(0) : k»> 



( ~ l 



\2N, 



c 



X" 



+ (7.75) 5 a (0| : (^7 CT 75y^)(0) : K(p)> 



2N< 
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<0| : (D p ^)?(0)V^(0) : |tt c ( P )) 



(z 75 )/*<0| : p p ^ 75T ^(0) : |tt c (p)) 



2iV, 



+(lal 5 ) s a (0\ : ( J D^ 7 CT 75y^)(0) : |vr»> 
+(7B^), a <0| : (D^i l5 ^)(0) : |vr c (p)> 
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2N C 
6N C J 



(0|:«0)(^W(0):K c (p)) 
'A d \ 



E 



A 



(* 75 y(0| : (i/ivn-DrfXO) : |vr c (p)) 



2iV, 



+(7.7 5 )/(0| : (^ 7 <T 75y^)(0) : |tt c (p)> 



2iV, 



+(75<V)/<0| : (^75y^»(0) : K(p)> (^-L 

In the above, we have also made use of the equations of motion in the chiral limit, pip = 0- 
Invariance under space-time translations, parity and charge conjugation further yields 

(0|(^75y^)(0)|vr c (p)) = (0|(^ 75 y^)(0)|vr c ( P )) = % - Vp <MV C , 



(0\(D p ij lal5 ^)(0)\Tr c (p)} = (0|(^75y^)(0)|vr c (p)) = ^ PpPr7 F 5 dc . 



Using these expressions, we deduce from Eq. ( |A.5| ) the following short-distance behaviour 
of the vertex functions 



aba 
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lim (rW)^ 
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feVP" + (P " ^) 2 ^ - (p • ?)(?mP" + SWV)] 
2 (p-q)q^-q 2 p^ 



1 W>)o , afec 

A? 2 F 7 



Ag 2 

v 



9/x + 



Ag 2 

p-q \ 
V J 



+ 



Ag 2 



(A.6) 



The result for IVp contradicts the one given in Eq. (55) in Ref. |33] where no term of order 1/A 2 
appears |)7j]. Note that this term is transverse, in accordance with the chiral Ward identities 
from Eq. (|2.2|) . On the other hand, our result for the short- distance expansion of Tva agrees 
with the one given in Eq. (58) in Ref. ||33|| . 
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B The resonance Lagrangian 

The resonance Lagrangian from Ref. [j38 |, which generalizes the one already given in Ref. [20 
reads (omitting terms including a scalar resonance) 



r 

*— re 



C 



(2) 

vv 



C (2) 



AA 



,(2) 
-VA 



Ly + La + £-vv + + £^ 2> 



VA i 



+t6 v e IIuaf3 (Wu»u a u 13 } + h v e^(W {u v , ff}) , 

-\{A^ V - 2Mi4>) - -Lf A (A,jn + ta A (A, [u„, /f ]) 

+7? (A^uW) + 7 ? (4 R, u"u„}) + 7 i 3) (Apu v ) (u>V) + 7 i 4) (AX) (n v u v ) 
+h A e^ a0 (A»{u» ) f^}) ) 

+hP(vxv u u" + %u v v v u») + \k v (%v»x+) + \iM% [Vu, /r]> 



1 



1 



iA^(% [A v J^])+iA^(% [u v ,A^]) + iA®(A fi [u v ,V^]) + B{% 
+He tlvaP (W {A\ ff}) + i^e^,( U V V^}> + zZ< a > 



(B.l) 

where the vector fields describing the vector and axial vector resonances have been denoted by 
and Ap, respectively In Eq. ( |B.1| ) we employed the usual notations |T8|, |2"U 



1 8 • 
d^R+ T^,R 



R = V,A 



]- {u ] [dp - + a M )] u + u [9^ - - a M )] m 1 "} 



uFPu< ± w f F 



R U i 
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: % [dp - i(v^ + a M )] u - u [9 M - - a M )] n f | = iu ] DJJv) = u\ , 
X± = « t X« t ± uy^ u , 



where the symbol (M) stands for the trace of the 3 x 3-matrix M. The field u is the square 
root of the Goldstone boson field, U = u 2 , whereas v^, a M and x denote the external sources. 

Note that in the case of bilinear interactions of and A M , only terms with at most one trace 
in flavour space are included in C res , which is compatible with the large- N c suppression of Zweig 
rule violating contributions. Furthermore, no pseudoscalar resonances have been incorporated 
in the resonance Lagrangian in Ref . |38| . As stressed in Ref . pO] , terms from £ 2 + £ 4 involving 



Goldstone bosons have to be added to the resonance Lagrangian as well in order to correctly 
reproduce the QCD short-distance behaviour of certain Green's functions. 
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